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Abstract. By using a recently proposed probabilistic approach, we determine 
the exact ground state of a class of matrix Hamiltonian models characterized 
by the fact that in the thermodynamic limit the multiplicities of the potential 
values assumed by the system during its evolution are distributed according 
to a multinomial probability density. The class includes i) the uniformly fully 
connected models, namely a collection of states all connected with equal hopping 
coefficients and in the presence of a potential operator with arbitrary levels and 
degeneracies, and ii) the random potential systems, in which the hopping operator 
is generic and arbitrary potential levels are assigned randomly to the states with 
arbitrary probabilities. For this class of models we find a universal thermodynamic 
limit characterized only by the levels of the potential, rescaled by the ground-state 
energy of the system for zero potential, and by the corresponding degeneracies 
(probabilities). If the degeneracy (probability) of the lowest potential level tends 
to zero, the ground state of the system undergoes a quantum phase transition 
between a normal phase and a frozen phase with zero hopping energy. In the 
frozen phase the ground state condensates into the subspace spanned by the states 
of the system associated with the lowest potential level. 
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1. Introduction 

A multitude of evolution problems can be cast in the form of linear flows, dfip = Hip, 
where is a matrix operator not necessarily Hamiltonian. For an initial value ip{0), 
the solution ip{t) of these systems of linear differential equations requires, as it is well 
known, the evaluation of the exponential of the operator H. The solution, with a real 
or an imaginary time t, also admits an exact probabilistic representation in terms of 
a proper collection of independent Poisson processes P El E] ■ 

Recently, we have exploited this probabilistic representation to derive exact finite- 
time solutions of a Hamiltonian flow by a Monte Carlo algorithm 0], as well as 
analytical results in the long-time limit OE]. In the latter approach, the Hamiltonian 
operator is decomposed, in a chosen representation, into diagonal and off diagonal 
parts, the potential and the hopping terms, respectively. On applying a central 
limit theorem to the rescaled multiplicities of the values assumed by the potential 
and hopping terms in the configurations visited by the system during its evolution, 
we have obtained a simple scalar equation whose solution provides an approximated 
semi-analytical expression for the lowest eigenvalue of H, in the following often called 
ground-state energy. 

The scalar equation derived in [2| contains only the first two statistical moments 
of the potential and hopping multiplicities and suggests that it is a second order 
truncation of an exact cumulant expansion. We have obtained this exact cumulant 
expansion by a large deviation analysis of the relevant probability density 0. In 
principle, were all the cumulants known, we would be in possession of a scalar 
equation whose straightforward solution is the exact lowest eigenvalue of H. In 
general, however, only a finite number of cumulants is available so that we have 
an approximated truncated equation whose level of accuracy depends on the system 
considered. In some cases, the convergence as a function of the number of cumulants 
is rather fast and the use of the first few cumulants provides results indistinguishable 
from those obtained by exact numerical simulations, see |Zj for details. 

In this paper we revisit the probabilistic approach El from a different point 
of view. We suppose that the asymptotic probability density of the potential and 
hopping multiplicities is known and derive an exact equation for the ground state 
in a closed form, i.e. not as a series expansion in cumulants. This formal result is 
then made concrete by observing that there exists a class of systems whose associated 
probability density is a multinomial. In fact, in the case of uniformly fully connected 
models and for random potential systems we obtain a very simple equation which 
relates the lowest eigenvalue of H, Eq, to that of the Hamiltonian operator with zero 
potential, Eq^\ This equation becomes exact in the limit M — > oo, where M is the size 
of the matrix representing H, provided that in the same limit E'q"'' diverges, namely 
a thermodynamic limit. 

By reason of the analytical expression found for Eq{M) as a function of Eq^\m), 
supposed known, the thermodynamic limit of the ground state of the above systems 
can be exactly characterized. We realize that, in this limit, a singularity may show up 
in the solution for Eq and a quantum phase transition takes place. Such a behavior is 
obtained when the degeneracy, or probability, the name depends on the kind of system 
we consider, of the lowest potential level vanishes. We also provide the equation 
of the critical surface separating the two phases and show that the hopping energy 
represents an effective order parameter of the transition. For both the uniformly fully 
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connected models and the random potential systems, the thermodynamic features of 
the ground state, including the existence of a quantum phase transition, are universal, 
in the sense that these depend, up to a rcscaling by |i?o°''li o^^y the levels and 
on the corresponding degeneracies (probabilities) of the potential. In particular, for 
the random potential systems we obtain a universal behavior independently of the 
nature of the hopping operator. This conclusion compares rather well with the results 
of exact numerical solutions for finite-size systems of quantum particles in one- and 
two-dimensional lattices. In fact, independently of the range chosen for the hopping 
operator, long or first neighbor, and independently of the nature of the particles, hard- 
core bosons or fermions, in all cases we have a clear tendency toward the universal 
thermodynamic behavior predicted by our formula. 

The paper is organized as follows. In Section El we introduce the uniformly fully 
connected models and the systems with random potential we deal with in the rest 
of the paper. In Section |31 we summarize our main result, namely a scalar equation 
whose solution provides the exact ground-state energy of the Hamiltonian operator 
for the above models in the limit M ^ oo. The proof of this result is given in Section 
0] In Subsection 14.11 we revisit the probabilistic approach in the case in which the 
asymptotic probability density of the potential and hopping multiplicities is known. 
The lowest eigenvalue of H is obtained, together with the asymptotic frequencies 
associated with the potential and hopping multiplicities, as the solution of a general 
system of equations. This result is specialized to the case of uniformly fully connected 
models and to random potential systems in Subsections 14.21 and 14.31 respectively. In 
Section [3 we discuss the thermodynamic limit and show its universal behavior and 
the appearing of a quantum phase transition when the degeneracy (probability) of the 
lowest potential level vanishes. Comparisons with exact numerical results obtained for 
finite-size lattice systems are provided in Section El in the case of random potentials 
with discrete spectrum and in Section [3 for a potential with continuous spectrum. 
The general features and results of our approach are summarized in Section (HI 



2. Uniformly fully connected models and random potential systems 

Let us consider a finite set of M states labeled by a vector index n. We will indicate 
this set as F, i.e. ¥ = {n} and |F| = M. The nature of the vector n depends on 
the context; in the case of an Hamiltonian particle model, for example, n represents 
an element of a Fock space. Let us also consider a state function V : ¥ ^ R with 
V{n) = Vn. We look for the lowest eigenvalue, Eq, of the M x M matrix H whose 
off-diagonal matrix elements are all equal to —77, with 77 > arbitrary, whereas the 
diagonal terms are given by the values Vn, with n e F. In equations, 

H = K + V, (1) 

where the elements of the hopping and potential matrices, K and V , are given by 

Kn.n' = -?7(1 - S.n.,n'), (2) 
Vn.n' ^ Vn^Yi,!!,' • (3) 

In the following, we will indicate with Y the set of all the possible different values of 
Vn with n G F, the levels of the potential for brevity, and with pv the degeneracy of 
level V, 
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Note that \Y\ < M and Y.vg'tPv = 1- 

The matrix H is the F representation of a Hamihonian operator H describing 
a fully connected model, or complete graph. This can be understood as follows. In 
general, a linear operator H defines the time evolution of the state function ip{n;t) 
of the system according to the equation dtipin^t) = Hip{n;t) starting from some 
given initial condition '(/;(n;0). We can always split H in two terms, K and V, such 
that in the F representation the corresponding matrices have, respectively, only off- 
diagonal and diagonal non vanishing elements. As it is clear from the probabilistic 
representation of the evolution equation, see Pj , the off-diagonal terms of H arc the 
rates for the transitions n n' between two different states, whereas the diagonal 
terms represent weights associated with the permanence in the states of the system. 
Therefore, Eqs. (|II3|I describe a system in which all the possible M(A/ — 1) transitions 
n — > n', with n ^ n' , take place with the same rate, namely a uniformly fully 
connected model. 

The other class of models which we will consider in this paper are defined in the 
following way. Given an arbitrary hopping matrix with elements Kn^n'-, the values of 
the potential matrix Vn are assigned to the states n randomly with the probability 

P(K. = T/)4^^' ''''/^ (5) 

where, now. Y is the set of the chosen levels of the potential V and {pv} the set of 
the corresponding given normalized probabilities, 



E 



PV = 1- (6) 

The choice of Y and {pv} is arbitrary. 

The two classes of models described above are, in a sense, complementary to each 
other. In fact, for the former we have the particular hopping operator (j^), which 
allows uniformly random connections among all the states, and an arbitrary potential 
operator, whereas for the latter we have an arbitrary hopping operator and random 
potential levels. 

3. Main result 

Consider first the uniformly random connected models. For V = 0, the lowest 
eigenvalue of which we call Eq^\ is trivially related to the lowest eigenvalue of the 
unit matrix and is given by 

4") = _^(A/-1). (7) 

It coincides, up to the factor —77, with the number of states connected to any state 
of the system. Similarly, for V constant, let say Vn = Vq for any n £ F, we have 
Eq = —rj{M — 1) + Vq. Note that, for the particular value Vb = g^t the obvious 

result Eq = —rjM. 

In the case of an arbitrary potential operator V , an expression for Eq is not 
known. By using the probabilistic approach developed in El j in this paper we 
will show that; in the limit A/ cx), Eq is the unique solution of the following equation 

Pv 1 



E 



Eq < Knin, (8) 
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or, in terms of the noninteracting energy iJg"'' , 

E ^ = ^0 < K... (9) 

It is easy to check that Eq. Q with Eq as an unknown has always solution and 
that the condition Eq < Vmin, where Vmin is the smallest element of ensures the 
uniqueness of this solution. 

Similarly, we will show that for M oo the ground state of the random potential 
systems defined by Eqs. (01 and © is also given by Eq. lO, provided that in the same 
limit E^°^ diverges. For these systems, the noninteracting ground-state energy has no 
more a trivial form like in ijJJ), as it is the lowest eigenvalue of an arbitrary hopping 
matrix K. In this case, E'g"'' must be determined as a separated problem, e.g. by using 
Monte Carlo simulations in the absence of sign problem, or by Fourier transformation 
if K represents the hopping operator of a system of independent particles. 

It is worth to stress that even if, in general, Eq is given only implicitly by Eq. 
(jHJ, this is a simple one-dimensional equation which can be straightforwardly solved 
by an iterative method. 

As we will show in details in Section 5, in the thermodynamic limit the solution of 
Eq. ^ may develop a singularity in its derivatives with respect to the potential levels 
y, signaling a two-state quantum phase transition with a frozen phase having zero 
hopping energy. These results hold for potentials with both discrete and continuous 
spectrum. 



4. Proof of Equation © 

The derivation of Eq. (jSJ follows from an analytical probabilistic approach we have 
recently developed to analyze the ground state of an arbitrary Hamiltonian operator 
H represented by a finite matrix. In subsection l4. II we review the basic definitions and 
results of the above approach referring the reader to [7| for details. Although in [7] we 
have explicitly considered lattice systems, i.e. the states n are lattice configurations 
with n indicating the occupation numbers of the sites, nothing changes if one looks 
at the states n as arbitrary abstract states. Moreover, differently from we will 
assume that the probability density which is the core of our approach is known. As a 
result, we will find a formal equation for Eq in a closed form. 

In general, given an arbitrary Hamiltonian matrix H, each row n of the 
corresponding hopping matrix K has a different number of non zero elements, let 
say A{n). We call A{n) the number of active links of the state n. When the number 
of active links of each state are all equal to the maximum value allowed, namely 
A{n) = M — 1, we recover a fully connected model. First, we will consider the 
analytical probabilistic approach in the general case. Applications to the uniformly 
fully connected models and to the random potential systems will be discussed in 
subsections 14 . 21 and 14.31 respectively. 



4.1. Analytical probabilistic approach 

Given the Hamiltonian operator H and separated its corresponding F-representation 
matrix H into the hopping and potential matrices, K and V, respectively, we define 
a virtual dynamics as follows. Let us parametrize the matrix K as 



Exact ground state for a class of matrix Hamiltonian models 



6 



such that |A„_„'| can be either or 1 and rjn^n' > 0. In graph theory the matrix with 
elements lAn^n'l is known as the adjacency matrix: in fact, it estabhshes whether two 
given states n, n' are first neighbors or not. We consider the Markov chain defined by 
the transition matrix P with elements 



A{n) 



(11) 



where 



A{n) ^Y.\^r^,r^'\■ (12) 



Starting from a given initial configuration no, we draw a new configuration rii with 
probability Pna.m- By iterating this procedure for TV steps we construct a path, or 
trajectory, in the space F n-o, rii, . . . , njv- For simplicity, in this paper we consider 
iJ to be a Hamiltonian operator, i.e. H \& a, complex Hermitian or real symmetric 
matrix. The approach can be generalized to non Hamiltonian operators. 

We will show that the information about the ground state of H is contained in the 
ensemble of the infinitely long paths. Along each finite path with N steps we have the 
sequences of data Aq, Ai, . . . , A^, Vq,Vi, . . . , Vat, Ai, . . . , Aat, and 771, . . . , Tyjv, where 

Ak = A{nk), k = 0,...,N, (13) 

Vk=^V{nk), fc = 0,...,7V, (14) 

Afc = Ajij._ J , k^l,...,N, (15) 

Vk = Vn.,^i,n,, k = l,...,N. (16) 

For later use we also define the sequence of values Tq, Ti, . . . , T/v-i, where 

Tk = AkVk+i/e, k = 0,...,N~l, (17) 

e being an arbitrary reference constant, which has the same dimensions of the yy's, 
typically an energy. Note that by construction |Afc| = 1. Let us indicate with Y, ^ 
and ^ the sets of all the possible different values that can be taken by the functions 
V, T and A, respectively. Let my, and be the cardinalities of these sets. In 
Uj we have shown that the ground state energy of H can be expressed in terms, not 
of the detailed sequences of the functions V, T and A, but just of the ensemble of their 
multiplicities, i.e. the number of times, iVy, Nt and N\, a given value for V , T and 
A has, respectively, occurred along each path, i.e. 

N 

iVy = ^<5y,y,, (18) 
fe=0 
N-1 

Nt^J2^t,t„ (19) 

fe=0 

JV-1 

iVA- 5]<5a,a.. (20) 

k=0 

Note that these multiplicities arc normalized to the number of steps N 

N+1, y = y - iv. (21) 
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More precisely, in |Zj we have proven that the matrix elements of the evolution operator 
at time t have the following probabilistic representation 

J2{n\e-"'\no)^Y.(wNi{Nv};t)l[T^-l[X^^\ , (22) 

n N=0 \ T<£S X<£^ I 

where (•)^ means an average over the paths of length N generated by the Markov 
chain starting from the initial configuration no. These averages over paths of 
fixed length are named canonical. The path functional weights Wat that appear in 
Eq. (|22|l are defined by the system 

WN{{Nvy,t) = / .2^, — ' (23a) 

^— = t, Xo > -Knin. (23&) 

For brevity, in the following we will drop the dependence of the weights on the potential 
multiplicities. 

Once we have access to the canonical averages and the series in the r.h.s. of Eq. 
H22|l is summed, wc obtain the lowest eigenvalue of H as 

Eo = lim -atlogV(n|e-^*|no), (24) 

t~^oo ^ ^ 

n 

with tlq arbitrary provided that it has a non zero projection onto the ground state 
\Eo). 

To evaluate the canonical averages it is useful to introduce the frequencies 
uy = Nv/N, V & vj, = Nt/N, T & .7, and = Nx/N, X e ^, which for 
N large become continuously distributed in the range [0, 1] with the constraints 

E '^^'^ = E = E = 1- (25) 

ver Tesr xe^ 
Note that, for N large, we do not distinguish the different normalizations, + 1 for 
Ny and N for Nt and N\, respectively. When possible, for the multiplicities and the 
frequencies we will use a compact notation in terms of the vectors fi and u, which have 
m = my + + components indicated by a Greek index a e = U U 
and are defined as = {... i/y Vj. i/-)^ ...) . We have 

fj, = Ni^. (26) 

For later use, we also define it"^ = (. . . — log[(a;o + V)/e] . log T . log A . . .), 
v'^ = (...(a;o+T/)-i...;...0...;...0...)andu;'r = .{xa+V)-^ .0 .0 .. .). 

Note that the vectors u, v and w depend on v through xq = xo{i') and v = —dx„u, 
w ~ —dxoV. Finally, we will take advantage of a scalar product notation. For instance, 
we rewrite Eq. (123 which determines xq, as {u, v) = t/N. 

By using the above compact notation, we express the iV-th term of the series in 
the r.h.s. of Eq. I|22(l in the following explicit form 



N 

where the probability VNifJ-) is given by the fraction of trajectories branching from 
the initial configuration Uq and having, after N steps, multiplicities fi. 
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In the hmit of long times t the paths most contributing to the evohition of the 
system are those with N ^ t large, see [Z] for details. In this limit, therefore, we can 
change from discrete paths with multiplicity fi to continuous paths with frequency u 

I \ r xot+N{u.u) 

(W^it) n T''^ n ^""O = / diNu)r^{N,.) 

N[xat/N+(v,u)+N-^ XogVN^Nu)] 

■\/27rA'e"^(i>', w) 

Here, we assume that the constraints (|25|) arc automatically taken into account by 
the probability VNiNv). Later, we will find more convenient to relax this feature 
of Vn{Nv) and explicitly introduce proper Lagrange multipliers. For t and N large, 
with N ^ t, we can evaluate the integrals in Eq. ()28(l by steepest descent 



yVN{t)Y[ T'^^ Yl X^"-) ~ e^«''*+^[('''''^"''')+^"''°s^"(^'^''')], (29) 

where u^^ is the saddle point of the non-smooth, exponentially-varying part of the 
integrand in Eq. H28|l and we added the superscript ''p to any function of u to indicate 
the value of this function for v = v^'p . The symbol ~ means asymptotic logarithm 
equality. The saddle-point frequency is the solution of the system of equations 

uaixo^u)) + N-^d^^ \ogVN{Nu) = 0, a e Jf", (30) 

which are obtained by derivating the exponent in Eq. (|28|) with respect to Va and 
using the property [u^v) = t/N. 

The series of the canonical averages in Eq. (|22|l is easily summed by replacing 
the sum with an integral over which is asymptotically exact for t — > cxd, see 0, 
and computing the integral by steepest descent. The result coincides, in the sense 
of asymptotic logarithm equality, with the r.h.s. of Eq. (|29|) evaluated at the saddle 
point N^^{t) solution of the equation 

■° [t- Ar(iy''P,t;^P)] 



dN 

= 0. (31) 



dN 

aeH Ar=7V=P(t) 

The last two terms of the above equation take into account the non explicit A^- 
dcpendcncc of the exponent in Eq. H29(l through XgP and u^^, respectively. We get rid 
of the third term in Eq. (|3H) by observing that for large we have 

lim N~HogVN{Nv) ^ uj{v), (32) 

where uj{v>) depends only on the frequencies ly, sec [Appendix A| The fourth and 
fifth terms in Eq. H31(l also vanish due to the property (i^'^p, d'^p) = t/N and the 
saddle-point equation for i/'^p, respectively. In conclusion, N^P{t) is determined by 

iu^P,u^P)+N-HogVNiNu^n\M=N''.it) = 0- (33) 
As expected from a physical point of view, wc have N^P{t) ^ t when t — > oo, see [7| 
for details. Finally, we obtain 

5^(n|e-^*|no) ^ e<^* (34) 
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Observing that for t large we must have 
energy Eq is given by 



-Hti 



no) 



hm 

t — >oo 



^0 \N=N^P{t) 



the ground-state 



(35) 



Taking into account the definition of u, which depends on xq, in the hmit t — > oo Eq. 
H33|l can be read as a a time independent equation that determines the quantity Eq. 
Note that in the same hmit i>'''^|^r=jvBp(() assumes a constant value. 

The key ingredient of the approach outlined above is the knowledge of the 
probability 7'7v(m)- Once this is known, the ground state can be derived analytically 
by solving the equation 

'-'\ogVN{N,.^n]j,=j,..^t)=0. (36) 



where lim* 



t 



lim I 

t — *oo 



JV=Af=p(t) 



is a finite vector dependent on Eq and 
(. . . - log[(-So + V)/e] . .logT. . . ; . . .logA. . 



Ar=Ar=p(t) 

T 



(37) 



Equation (|36(l allows the determination of the ground-state energy Eq in terms 
of asymptotic saddle-point frequencies, for brevity v hereafter, which are determined 
by the system (|30|l in the limit N ^ oo. Therefore, by using again Eq. H32(l . we find 
that Eq and v are the solution of the following system of coupled equations 



{v>,u{Eq))+uj{v) 
H ^v- = 1, 

E 
E 



0, 



0, a e jr, 

Eq < Vmin, 



(38) 



1, 



1, 



where v?- = (. . 



log[(--Bo 
cr, 



< 



csr, 



V)/e]...-. 

a e r. 



■logT. 



. log A . . .) and 



(39) 



are the three Lagrange multipliers introduced to explicitly take into account the 
constraints (|25|l so that the function ijj(i'), i.e. the corresponding probability, is now 
meant to have support on the whole hypercube [0, 1]™. In other words, the function 
wlu) in the system (|38|l is the analytic continuation to the set [0, 1]™ of the function 
defined in Eq. 

In |2] we have used a cumulant expansion theorem to express, in the general case 
in which VnIp) is not known, Eq as a function of the cumulants of the variables V, T 
and A generated the Markov chain over the states n. It is clear that in this case only 
the first few cumulants can be considered available, numerically or analytically, and 
the expression found for Eq is necessarily approximated. On the other hand, the exact 
value of Eq can be obtained in the cases in which VnIp) is known. As we shall show 
in the next two subsections, the uniformly fully connected models and the random 
potential systems provide non trivial and important examples of the latter class. 
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4-2. Uniformly fully connected models 

Up to now we have formulated the problem in the most general case. In this Section, 
we specialize the discussion to the uniformly fully connected models defined by Eqs. 
H1I3|) . For these models, both the sets ^ and Jif have a single element, namely 
T = {M — l)'r]/e and A = 1, whereas we may count, in general, M distinct values V in 
the set y. This implies that we have to consider only the potential multiplicities 
Nv, V e i.e. (1^ = {. . . Nv . . .■,N;N). Furthermore, due to the fact that 
all the transitions n — > n', with n ^ n' are allowed with the same probability, 
we can assume that after N steps the Markov chain Hll() is well represented by a 
multinomial distribution, in which, at each step, any state has the probability 1/M 
to be extracted. This representation will become exact in the limit M oo, i.e. 
when the introduced extra transitions n ^ n become immaterial. According to 
this multinomial distribution, the probability to have, after N steps, multiplicities /x 
satisfying the constraint X^yey = is given by 

^-('^) = ^!n^- (40) 

ver 

The parameters pv, with V ^ y, are the degeneracies of the potential values V in 
the space of the states F, see Eq. Q. Due to the ergodic properties of the Markov 
chain pi(l . they also coincide with the asymptotically expected frequencies of the 
multiplicities Ny 

p. = ^lm^l(iV,>,. (41) 

The latter expression is relevant to calculate by Monte Carlo simulations the 
parameters pv in those cases in which a direct evaluation in the space F is not feasible. 

By using the Stirling approximation for N large, the probability for the potential 
multiplicities iVf, with each component of v continuously distributed in the range 
[0, 1], can be written as 



Vn{Nu) ~ exp 



V vv log ( — ) 



N 



(42) 



from which we get 



u{v)^Y.^y^^^{—)- (43) 

Taking into account that for this model we have VT={M-i)ri/e = 1 and h'x=i ~ 1, 
the only random multiplicities being those of the potential values, the system of Eqs. 
H38|l becomes 

uv(i;o) + log(^^^ -1 + cr = 0, 

^ vvuv{E,,) + log [(A/ - l)77/e] + log (^-^ =0, < Kni^, (44) 
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By using uv{Eo) = — log[(— E'o + V)/e] and observing that the Lagrange multipher 
cy does not depend on the index V, from the first equation of the above system we 
find 

where the normalization constant, Z = e^^cxp(c:^ — 1), is determined by using the 
third equation in H44|) . namely 

By inserting the result (|45|l into the second equation of the system (I44|) . we get 

1 log(Ze) + log [(M - l)Ve] = 0, (47) 

which, on using the value H46(l for Z, brings immediately to Eq. (jHJ. 



4-3. Random potential systems 

The above derivation can be readily extended to those systems in which, both the 
hopping values T and the phases A are general, whereas the potential levels V € Y 
arc independent random variables assigned to the states n with arbitrary probabilities 
Pv, sec Section 121 It is clear that different realizations of the values Vn correspond 
to different Hamiltonian matrices, i.e. we deal with an ensemble of random matrices. 
In this case, therefore, the study of the ground-state energy should be meant as the 
evaluation of the average lowest eigenvalue of the ensemble and of its fluctuations 
around this value. In the following, we will assume that the probability for the 
ground-state energy to assume a given value becomes infinitely peaked around its 
average Eq. 

For a finite system with a specific realization of the potential values Vn the 
corresponding ground-state energy is affected by the correlations among the Vn^s 
induced by the virtual hopping dynamics. Consider, in particular, the possibility 
that a state n comes back to itself after a few steps. However, such correlations must 
vanish if the number of active links associated with each state, A{n), diverges. In 
this limit, any step of the Markov chain is equivalent to a random extraction of the 
value V distributed according to the given py. In fact, the situation in which the 
numbers A(ri) increase by increasing the number of states M is rather common, later 
we will consider the important example of many-body systems on a lattice. For these 
systems, in the limit AI ^ oo the potential becomes exactly uncorrelated from the 
states n and along the virtual dynamics there holds the following factorization 

Pn{{Nv},{Nt}ANx})^P^{{Nv}) V^^\{Nt}ANx}), (48) 

where the superscripts ^ and '^"^ refer to the probability density for the potential 
multiplicities {Ny} and for the rest of the variables {Nt} and {A^a}, respectively. In 
terms of the asymptotic function uj defined by Eq. (|32() , this implies that 

^(«y) = u^iM) + uj^'HM, {^a}). (49) 

The consequences of the factorization (|48|) are immediately obtained. Note that in 
the probabilistic representation illustrated in Subscction l4.1l unlike other probabilistic 
representations in which importance sampling is included, the frequencies vt and vx 
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are independent of the potential values, so that their distribution can be evaluated for 
V = 0. In this case, we have i'v=a = 1 and uv=o ~ ^ log(— ^^q^V^) and, therefore, 

)) = -log(-4"Ve) + E ^TlogT+ E ^.AlogA. (50) 

From this relation and from the second equation of the system H38|l . we find that the 
energy of the nonintcracting system, Eq^\ satisfies the equation 

log(-4"Ve) = E '^Tlogr+ ^^AlogA + t^W({i.T},{^A}). (51) 

From Eq. 1)49(1 and assuming that E^^^ is known, we can use Eq. H51|l to get rid of 
the quantities in the r.h.s. in the general case V 0. In fact, by inserting Eqs. ((49|1 
and 1(5 1() in the general system 1(38(1 . we are left with the following reduced system 

(uviEo) + d,^uj'^' {{vv]) + cr = 0, 1/ S r, 
E vvuv(Eq) + log(-i;<"Ve) + ^'^{{vv}) - 0, £^0 < K.i„, 



(52) 



Finally, by using the fact that the potential levels F G 1^ are randomly assigned with 
probabilities py, we have that uj^ {\yy\) is again given by the multinomial formula 
1(43(1 . Therefore, apart from the exchange of E^^ with — (M — 1)77, the system ((52(1 is 
formally identical to the system ((44(1 so that, in the limit M 00, the ground-state 
energy Eq of the random potential systems satisfies Eq. . 

5. Thermodynamic limit: quantum phase transition and universality 

For the uniformly fully connected models and the random potential systems, the 
solution Eo{M) of Eq. approaches the exact ground-state energy when the 
number of states M becomes infinitely large provided that the number of active links 
associated with each state, A{n), also diverges. The condition on the number of active 
links, which is evidently fulfilled by any uniformly fully connected model, for random 
potential systems can be viewed as a condition on the noninteracting ground-state 
energy Eq^\m), which is proportional to a proper average of A{n) with n G F. In 
this Section, we will study a non trivial thermodynamic M 00 limit of Eq. @ 
always assuming that lim^/^oo \Eq^\m)\ ~ 00, as happens for any system of particles 
with fixed density. 

For simplicity, we suppose that the number of potential levels, i.e. the 
cardinality my of the set Y, does not change by increasing the number of states 
M. Consistently, we assume that the degeneracies (probabilities) {pv{M)} tend 
to constant values for M — s- 00. Looking at Eq. lO, it is immediately clear 
that if the potential levels {V{M)} diverge for Af — > 00 more slowly than 
E^°\m), we have limj^f _oo -^o (^)/^^o°^ (^) = 1- Another trivial hmiting solution 
is found if the potential levels {V{M)} diverge faster than Eq^\m), namely 
lirRM—Kxy Eq(M) /VminiM) = 1. A thermodynamic non trivial limit is obtained only 
for V{M) ^ E^q\m). In this hmit, we rewrite Eq. I© as 

^" =1, eo < Wmin, (53) 

^ - eo 
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where we have defined py = Iuhm^oo Pv{m), v = liinM^ooV{M)/\EQ^\AI)\ and 
Co = limM^oc Eq{M)/\Eq^\m)\ and we have assumed, as usual if K represents a 
hopping operator, that Eq^\m) < 0. With an abuse of notation, we keep using 
the same symbol chosen for the set {V} also for the set of the asymptotic rescaled 
potential levels {v}. The asymptotic rescaled ground-state energy eo cannot exceed 
the minimum asymptotic rescaled potential level Wmin = min„gt'{w}- In the following, 
we will find handy ordering the elements oi "f as vi < V2 < ■ ■ ■ < I'm^ , and indicate 
with pi,p2, ■ . ■ ,Pmy the corresponding degeneracies (probabilities). 

For what concerns the uniformly fully connected models, although Eq. (|53|1 
is formally equivalent to Eq. it has a reacher structure. In fact, unlike the 

finite M case in which the py's vary in the rational field, sec Eq. Q), in the limit 
M oo the Pv 's vary in the continuum so that we are allowed to consider the analytic 
continuation of the solution eg toward limits in which one or more of the py 's tend to 
zero. Regarding the random potential systems a similar observation, strictly speaking, 
does not apply. In this case, in fact, the py's can vary in the continuum even with a 
finite M, see Eq. On the other hand, as explained in the previous Section, for the 
random potential systems the ground-state energy becomes a well defined quantity 
with vanishing relative fluctuations only in the limit M ^ oo. 

We shall show that, due to the constraint eg < Wmin, for pi — > the solution eo 
of Eq. (|53|l becomes non analytic and a quantum phase transition takes place. 

First we illustrate in detail an example with only two potential levels. In this 
case, Eq. I)53|l is a quadratic equation for cq which can be solved explicitly. Observing 
that p2 = 1 — Pi , we find 



eo = wi - ^ 



V(«2 -Vi- 1)2 + 4pi{v2 - Vi) - {V2 - Vi - 1) , (54) 



the other solution being incompatible with the condition eo < ui. The behavior of 
eo — Vi as a function of pi and V2 — wi is shown in Fig. ^ For pi — > we have 



lim eo - wi = 5 [('^2 - wi - 1) - \v2 - I'l - 1|] 

_ j V2 - Vi - 1, V2-Vi< 1, 

I 0, 'i;2 - I'l > 1, 



(55) 



i.e. a singularity shows up at ?J2 — = 1, sec Fig. |21 It is trivial to check that for 
Pi = 0, Co has a discontinuity in the first derivative and a divergence in the second 
derivative with respect to t)2 — i^i at i;2 — t^i = 1. In conclusion, for pi = we have 
a quantum phase transition at the critical point V2 — vi = 1. We stress that this 
behavior emerges due to the constraint eo < vi, which keeps to hold only if the limits 
performed are taken in the correct order: M ^ oo first and pi ~> later. 

It is easy to generalize the above result to a case with rui/ > 2. Let us rewrite 
Eq. ^ as 

-^^+E7 r^, V = eo-v,<0. (56) 

The unknown eo — wi is a function of the 2(rny — 1) parameters Pk,Vk — wi, 
k = 2, . . . ,my. Note that pi is fixed by the normalization condition. We show 
the graphical construction of the solution of Eq. (|56|l in Fig. O For pi > 0, i.e. 
SlSPfc < 1, we always have eo < vi and, furthermore, eo is an analytic fmiction of 
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Figure 1. Thermodynamic limit of a system with two potential levels: 
asymptotic rcscalcd ground-state energy eo — vi as a, function of pi and V2 — vi, 
see Eq. ICT . 




0.5 1 1.5 2 

Figure 2. As in Fig. Q sections of eo — di as a function of V2 ~ vi for different 
values of pi. For pi = a quantum phase transition takes place at the critical 
point V2 — vi = 1. 



its arguments. When pi = 0, a singularity of Cq shows up, which can be characterized 
in the following way. Let us define the function 

m-i/ 

W{p2,V2 -Vi,... ,Pm-^,Vm^ - = ^ . (57) 

t^2 - ^1 
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From Eq. (|56|l we see that, if pi = and eo < wi, then 



W>yi e^^T r = l. (58) 

On the other hand, if eo — > ui for pi —i- 0, from the same Eq. H56|l we have 

l-W = hm > 0, (59) 

pi-»0,eo^Uj" Vl — eo 

i.e. W < 1. In conclusion, for pi ^ the below relations hold between W and eo 

VF>l^eo<wi, (60) 
< 1 4=^ eo = wi. (61) 

These two equations establish the following scenario for pi 0. As we move inside 
the 2{m-r — 1) dimensional region determined by the condition W < 1, the rescaled 
energy eo stalls at its maximum value eo = wi. Outside this region, we have eo < wi. 
Therefore, we find that Veo|j^+ ^ and Veol^^- = 0, where ^+ and ,9'~ are generic 
points arbitrary close to the surface determined by the condition W ~ 1, and such 
that W > 1 01 W < 1, respectively. More precisely, we find that on the critical surface 
=5^ any directional derivative of eo has a discontinuity for any direction not tangent to 
c5^, which, in turn, implies a divergence of any double derivative of eo for any direction 
not tangent to o5^. 

We observe that according to Eq. H57|) the critical surface is determined by a sum 
over all the potential levels of the ratios Pk/{vk — vi). Therefore, even levels very far 
from the lowest one, vi, can contribute to the critical behavior if the corresponding 
degeneracies are large enough. This is a cooperative phenomenon due to the intrinsic 
quantum nature of the considered systems. 

We wish to discuss now in more detail the nature of the phase transition 
corresponding to the above singularity and also determine an order parameter. Let 
us write the potential operator V in terms of the projectors onto the configuration 
subspaces at fixed potential values 

V£-r neF:V{n)=V 

On derivating the ground-state energy Eq with respect to a potential level V and using 
the Hellman-Feynman theorem, we find 

dEp ^ d {Eo\H\Eo) ^ {Eo\dvH\Eo) ^ 

dV dV {Eo\Eo) {Eo\Eo) ^ ' 

where 

{Eo\nv\Eo} ^ \{n\Eo)\^ 

{Eo\Eo) ^ ^ T^;por ^ ^ 

is the ground state expectation of the projector Try. On the other hand, in the limit 
M — ^ oo we have dyEo — d^eQ so that by derivating Eq. H53|l with respect to a generic 
V, we get, in the notation in which the potential levels are ordered, 

Oeo Pk 



Ck = 



dvk {vk - eo)^ 



7n f 

Pk' 



1, . . . , my. (65) 
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Figure 3. Graphical solution of Eq. 1561 in a system with four potential levels 
for pi = 0.1 (top) and pi = (bottom) in two cases W > I (left) and W < I 
(right). The continuous lines are the l.h.s. of Eq. 1561 plotted as a function of eg 
whereas the horizontal dashed lines represent the unit level. The rescaled ground- 
state energy is given by the unique intersection (dots) in the region eo < vi. For 
pi = 0, a singularity shows up at = 1, namely eg stalls at its maximum value 
eo = ^'1 for W < 1. 



Equation holds at any point of the 



Vi,...,P: 



Ci 

whereas for fc = 2 




2(mr — l)-dimensional space (p2,^2 
H6UI61|1 and (jHSJ we get 



Pk' 



W <l. 



W >l. 



(66) 



(67) 



{vk' - eoY 

We deduce that in the thermodynamic hmit we have the following behavior of 
the ground state \Eq). In the region W > 1, |£'o) turns out to be an analytic function 
of its arguments pk,Vk — vi, k = 2, . . . , my, while in the region < 1 it collapses 
into the subspace spanned by the configurations with minimum potential value Vi. 
The susceptibilities 9(„j^_,„j)Cfc diverge on the critical surface y determined by the 
condition W = 1. Finally, it is simple to check that the asymptotic rescaled hopping 
energy in the ground state, eo°P = \imM^cs,{EoiM)\K\Eo{M)) /\E^"\m)\, is given by 

-1 

(68) 



hop 



E 

.k=l 



Pk 



(wfc - eo)2 



Therefore, the phase for < 1 is a frozen phase with Bq"^ = 0, whereas a non 
vanishing hopping energy is obtained for W > I. In Fig. 0]we show the behavior of 
^^°^{Pij V2 — vi) in the case with two potential levels described above. For pi = 0, it 
is evident the discontinuity of ejj"^ at the critical point W = 1, i.e. W2 ~ ^'i = 1- 
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Figure 4. Thermodynamic limit of a system with two potential levels: 
asymptotic rescaled hopping energy eg°'' as a function of pi and V2 — vi. 



6. Many-body lattice models: random potential with discrete spectrum 

Hereinafter we will focus our analysis on the thermodynamic limit of random potential 
systems. Equation H53|) states that for these systems the asymptotic rescaled energy 
Co is universal, independently of the nature of the hopping operator, provided that 
hmM^oo l^^o"^ = oo. In this Section we present the results of numerical 

simulations which help to quantify the rapidity with which this universality is 
approached by systems with different complexity. 

To begin, we investigate quantum particles moving in two-dimensional lattices and 
interacting via a potential with two random levels. Three cases have been considered: 
spinless hard-core bosons with long range hopping, spinless hard-core bosons and 
spinless fermions with first neighbor hopping. The results are displayed in Figs. 
El and respectively. For these systems, in Table Q] we report the noninteracting 
ground-state energy Sq^-* and the average number of active links 



as a function of the dimension M of the Fock space. In the case of fermions, due to 
sign cancellations, the effective average number of active links which is responsible of 
the potential decorrelation is the difference 




(69) 




(70) 
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Figure 5. Universality of tlie thermodynamic limit of random potential systems: 
rcscaled ground-state energy eo as a function of the rescaled potential level V2 
for systems of spinloss hard-core bosons with long range hopping. The random 
potential has two levels: Vi = with probability pi = 0.6 and V2 = ^21^0"' I 
with probability p2 = 0.4. The solid line is the universal thermodynamic limit 
predicted by Eq. 1531 whereas the dots corresponds to the numerical results 
found for the finite-size systems indicated in the legend. Here, mxn.N connotes 
a system of N particles in a m X n lattice. The error bars represent the standard 
deviation of the stochastic variable eo as evaluated from an ensemble of 100 exact 
diagonalizations of the Hamiltonian matrix with the specified random potential. 



Table 1. For diff'erent systems mxn.N of N particles in a two-dimensional mxn 
lattice with periodic boundary conditions we report the number of states M, 
the noninteracting ground-state energy Eg'^^ and the average number of active 
links Aav in the case of spinless hard-core bosons with long-range (hcb.lrh) and 
first-neighbor hopping (hcb.fnh) and spinless fermions with first-neighbor hopping 
(f.fnh). 



system 


M 


hcb.lrh 


hcb.fnh 


f.fnh 














3x3.4 


126 


-20 


20 


-10.34699 


10 


-7 


3.33333 


3x4.4 


495 


-32 


32 


-12.10548 


11.63636 


-9 


4.84848 


3x4.6 


924 


-36 


36 


-13.71426 


13.09090 


-10 


4 


4x4.4 


1820 


-48 


48 


-13.25169 


12.8 


-10 


5.52087 


4x4.5 


4368 


-55 


55 


-15.30041 


14.66666 


-12 


5.64102 


4x5.5 


15504 


-75 


75 


-16.44602 


15.78947 


-13.23607 


7.22394 



between the active links A±{'n) associated with positive or negative jumps, 
respectively. 

The case of hard-core bosons with long range hopping presents the closest analogy 
with a uniformly fully connected model. The number of active links A{M), although 
being much smaller than the number of states M and diverging more slowly than 
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Figure 6. As in Fig. |S]in the case of systems of spinless hard-core bosons with 
first neighbor hopping. 




^2 

Figure 7. As in Fig. |^ in the case of systems of spinless fermions with first 
neighbor hopping. 



M for M — > oo, is a constant which docs not depend on the configuration n. In 
particular, we have that A{M) = A^^{M) = -E^q\m) is given by the simple formula 
(number of sites — number of particles) x number of particles. The numerical results, 
shown in Fig. Elby dots with error bars, have been obtained by exact diagonalizations 
of the Hamiltonian matrix for several realizations of the random potential. The error 
bars represent the standard deviation of the stochastic variable Eo/\Eq^^\ evaluated in 
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0.15 



0.1 



T" 



T" 



o hcb.lrh 

— - (-0.002+0.0 10)+(0.560+0.647)/A,^^ 

< hcb.fnh 

— . (-0.008+0.042)+(0.599+0.622)/A^^ 

o f.fnh 

- . . - (-0.01 9+0.027)+(0.509+0. 1 58)/A 



0.05 - 







0.05 0.1 



0.15 

1/A 



0.2 



0.25 



20 



0.3 



Figure 8. Difference Aeo between the rescaled ground-state energy ep predicted 
in the thermodynamic limit for V2 = 100 by the multinomial formula I53i and that 
obtained in the finite-size systems analyzed in Fig. |S] (hcb.lrh o), Fig. |H] (hcb.fnh 
<) and Fig. □(f.fnh o) as a function of the inverse of the average number of active 
links Aav The lines are a weighted linear fit with errors on Aeo assumed equal 
to the standard deviations of eo for the corresponding finite-size systems. 



this way. The agreement with the universal curve for eg, obtained from Eq. 1)53(1 as a 
function of the second rescaled potential level V2, worsen by increasing the value of W2- 
However, we have a clear tendency toward the universal result by choosing systems 
closer and closer to the thermodynamic limit. 

For systems with first neighbor hopping, we have a much greater complexity with 
respect to a uniformly fully connected model. The number of active links is a function 
of the configuration n and, in the case of fermions, a sign problem is also present. 
Therefore, we expect a convergence to the universal thermodynamic limit slower than 
that obtained with long range hopping. This is confirmed by the results shown in 
Figs. Eland[3 Whereas the trend toward the universal behavior is clear in both cases, 
for fermions the convergence is indubitably slower than for hard-core bosons. 

Note that for practical reasons, essentially the finite amount of memory (4 
Gb) of the computer used to perform the numerical diagonalizations, the different 
systems considered in Figs. |S| |H1 and |7| do not have exactly the same density, 
number of particles/number of lattice sites, as it would be convenient in looking at 
the thermodynamic limit of a lattice particle system. In fact, we have chosen a 
set of systems in which M and |£'o'^''|, i.e. Aav, possibly both increase compatibly 
with the condition M < Mmax, where Mmax is the size of the largest diagonalizable 
Hamiltonian matrix. At constant density this set would contain only two or three 
elements, so we decided to allow for density fluctuations. These fluctuations may 
reflect a non monotonous approaching to the thermodynamic limit. Such a behavior 
is quite evident in the fermion case of Fig. [7|where the system 3x4.6 with M = 924 
and = 4 is closer to the asymptotic universal curve than the next system 4x4.4 
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which has M = 1820 and A^^ = 5.52087. 

In Fig. IHl we show a quantitative test for the convergence of data in Figs. [3 El 
and 13 The difference Aeg between the rescaled ground-state energy of the finite- 
size systems and the thermodynamic value predicted by the multinomial formula (|53|l 
is reported, for V2 ~ 100, as a function of A~^. In all cases, the data are fitted by 
straight lines which, compatibly with the associated errors, give Aeo for — > 0. 
This behavior corresponds to the scenario depicted in Section 14.31 as the average 
number of active links diverges, the multinomial formula l|53|) becomes exact with 
a residual error proportional to A~^. This law also explains the progressively larger 
errors observed, at a fixed system size, passing from hard-core bosons with long range 
hopping to hard-core bosons with first neighbor hopping to fermions with first neighbor 
hopping, see Eqs. and f7n| . 

The models considered above, even if share some realistic features with systems 
of interest in physics, have been studied in connection with a toy random potential 
defined by only two levels Vi and V2 with assigned probabilities pi and p2- For these 
models, in order to obtain the phase transition in the ground state we must impose 
the additional condition — > after the thermodynamic limit has been taken. In 
the remaining part of this Section, wc discuss a more realistic random potential which 
has the following characteristics: i) the number of levels increases by increasing the 
number of states M, and ii) the probability associated with the lowest level vanishes 
for M — > 00. 

Consider a lattice A with some ordering of the lattice points 1,2,..., rt[A| occupied 
by N quantum particles. As above, we assume that the particles are fermions or hard- 
core bosons so that the occupation numbers can be n,; = 0, 1, i = 1,2,..., |A|. To 
each of the M = |A|!/(A^!(|A| - iV)!) Fock states rt = (ni, n2, . . . , n|A|) of the system 
we associate the potential 



where ri(ri) arc a set of |A|M independent random variables assuming the values 
0,1. Let p be the probability that ri(ri) = 1. Since the Fock states are normalized 
by the condition X]l=i = -^i have -I- 1 different potential levels = 
k — 0,1, . . . , N with associated probabilities 



Note that po, the probability of the lowest potential level Vq = 0, vanishes for N ^ oo 
as far as p > 0. 

Before going further, we want to spend a few comments about the introduced 
potential H71|l . It represents the on site random coupling of a set of |A| impurities with 
the particles of the system. The impurities, spins to be concrete, may be in one of 
two available states with probabilities p and I — p. Impurity- impurity interactions are 
neglected. Moreover, the potential 1)71(1 treats the impurities itself as classical spins, 
i.e. it does not take into account the correlations of the state of each single impurity 
with itself in the presence of different dispositions of the particles in the lattice (Fock 
states). The latter decoupling could be the effective result of an impurity dynamics 
fast on the time scale of the evolution of the particle system or of the presence of an 
external random field. 




(71) 



Pk = 



p\l~p) 



N-k 



(72) 



ki{N-ky. 
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Table 2. As in tabic ITl in the case of one-dimcnsional systems Ixn.N. 



system 


M 


hcb.lrh 


hcb.fnh 


f.fnh 














1x8.4 


70 


-16 


16 


-5.22625 


4.57142 


-4.82842 


3.42857 


1x10.5 


252 


-25 


25 


-6.47213 


5.55555 


-6.47213 


5.55555 


1x12.6 


924 


-36 


36 


-7.72741 


6.54545 


-7.46410 


5.45454 


1x14.7 


3432 


-49 


49 


-8.98792 


7.53846 


-8.98791 


7.53846 


1x16.8 


12870 


-64 


64 


-10.25166 


8.53333 


-10.05467 


7.46666 



In the thermodynamic hmit |A|,iV oo with iV/|A| eonstant, the rescaled 

ground-state energy cq = E()/\Eq^^\ is determined by Eq. H53|l . In this hmit the sum 

over the discrete rcscalcd levels Vk = Vfc/|i?g''''| can be approximated by an integral 
over k so that Eq. H53|) reads 

^ -dk = 1, eo < 0. (73) 



Vk - eo 

Since Eq. (|71|l is a linear combination of the independent random variables ri{n), in 
the thermodynamic limit the distribution of the rescaled potential values, w^, becomes 
infinitely peaked around the mean value v = wjvp, 

Pk ^ 6{k - Np), (74) 

which inserted into Eq. (|73|l gives 

^ 1, eo<0. (75) 



VNp - eo 

In conclusion, the solution for the rescaled ground-state energy is 

r -l + jNp/\E^"^\, lNp/\E^°^\ < 1, 
eo = < (76) 
lo, ^Np/\E^°^\>1. 

Whether or not the solution given above represents a non trivial rescaled ground- 
state energy depends on how fast the potential levels diverge in the thermodynamic 
limit with respect to the noninteracting ground-state energy. For a system of 
fermions with first neighbor hopping Eq^^ can be evaluated analytically by Fourier 
transformation. In one dimension with periodic boundary conditions, at half filling 
and neglecting spin, for — > cxo we have Eq^^ ~ —Air^^Nrj, where ry is the value of 
the hopping coefficient between any two first neighbors. Equation H76() then becomes 

[ -f + njp/4:T], njp < 477, 
U, TTjp > Ar]. 

A comparison of eo predicted by the above expression for p = 0.4 with the values 
obtained from numerical simulations of finite-size systems is shown in Fig. |5| The 
main parameters of the simulated systems are given in Table El The numerical 
results are consistent with a phase transition occurring in the thermodynamic limit at 
7/77 = 4/7rp~ 3.18. 

The value of the noninteracting ground-state energy can be given analytically also 
for spinless hard-core bosons with long range hopping. At half filling with periodic 
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0.01 0.1 1 10 100 

y/Ti 



Figure 9. Universality of the thermodynamic limit of random potential systems: 
rcscaled ground-state energy eg as a function of the rescaled potential strength 
■y/r] for systems, Ixn.N, of N spinless fermions with first neighbor hopping in a 
one-dimensional lattice with n sites. The potential is the impurity potential of Eq. 
1711 with p = 0.4. The solid line is the universal thermodynamic limit predicted 
by Eq. 1771 whereas the dots corresponds to the numerical results found for the 
finite-size systems indicated in the legend. The inset shows a linear fit of Aeo vs. 
as displayed in Fig. |S] 




Yj/T] 

Figure 10. As in Fig. Elin the case of spinless hard-core bosons with long range 
hopping. The potential is the impurity potential of Eq. II71I with p = 0.4 and 
7 = 71 Af, with 71 constant. 
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boundary conditions we have Eq^"^ = —N'^r], where r] is the value of the hopping 
coefficient between any two sites. In this case, for 7 constant the potential levels 
diverge for — > cxo slower than E'q"' and Eq. (|76|l always gives eo ~ ~1. A non trivial 
result is obtained assuming 7 = 71 A^, with 71 constant, 

[-I + 71P/7?, liP<V, 
eo = < „ . (78) 

[ 0, 7ip > 77. 

In Fig. ^1 we show the behavior of Eq. H78|) for p = 0.4 in comparison with numerical 
simulations for finite-size systems. Again, there is a consistent matching with the 
singularity displayed by eo at 71/77 = 1/p = 2.5. As discussed in Section |5l in the 
present case the approach to the thermodynamic behavior is faster (smaller differences 
Aeo) with respect to the case of spinless fermions with first neighbor hopping shown 
in Fig. El 

Note that, for clarity, we have plotted in Figs. 1^ and II Ul different curves for cq. 
However, as it is evident by Eqs. H77() and H78() . the behavior of the rescalcd ground- 
state energy is universal once it is expressed in terms of the rescalcd potential levels 
Vfe/|i?Q°^|, provided, of course, that the same probabilities pk arc used. 

7. Many-body lattice models: random potential with continuous 
spectrum 

The results obtained in Section are readily extended to the case of a random 
potential with continuous spectrum. In the thermodynamic limit, for a continuous 
distribution of rescaled levels v = lim^v/^oo ^ (M) / \ Eq^\m)\ described by the density 
p{v) = liTanj^oo p{y{M)), the equation determining the rescaled energy eo = 
limM^coSo(M)/|4°^(M)| reads 

-^^dv = 1, eo < Vn,in- (79) 

w - eo 

The only point which wc have to pay attention to concerns the definition of the lowest 
potential level. In order to avoid ambiguities, we define Umin as the value of v such 
that 

P(wmin + 5) >0 and ^(l^min - 5) = 0, (80) 

with (5 > arbitrarily small. Note that definition H80() does not include the case 
fmin = — oo, which occurs, for example, if the density p{v) is a Gaussian. However, 
since for Wmin = — oo we can only have the trivial result eo = I'min = — oo, in the 
following we will assume Wmin finite. 

In order to recover the analogous results of Section [S] it is useful to define 

Fix) = [ P^dv (81) 
J V - X 

and 

W = F(«^i„) = / -^^dv. (82) 

J V - Vnun 

Observe that F(x) is a positive monotonously increasing function for x G (— oo, Umin] 
with 

- oo < F(a;) < = sup x G (-00, fmin]- (83) 



Exact ground state for a class of matrix Hamiltonian models 



25 



We have that eo is solution of Eq. l(7(^ if F(eo) = 1 and cq < Wmin- If p{vmin) > 0, 
and therefore W =^ oo, Eq. (|83|l shows that the solution eo exists, it is unique and 
smooth, i.e. without singularities with respect to the parameters of the density If 
p(wmin) = then, since / p(v)dv = 1, we have W < oo. For 1 < < oo, we still have 
a unique smooth solution eo of Eq. (|79() , with eo < Umin for > 1 and eo = f min for 
W = 1. For W < 1, Eq. H79(l has no solution and, as discussed in the discrete case, we 
have to consider for eo its analytic continuation for p(wmin) 0. More explicitly, we 
first solve Eq. H79(l with a modified density p{v) having the same support of p{v) but 
with p{vniin) > and then evaluate eo as the limit for p — > p of the solution so found. 
Clearly, this analytic continuation cannot exceed Wmin and it is easy to demonstrate 
that we have eo = Wmin- In fact, if it were eo < Umin, there would be no singularity 
in the integrand of Eq. (|79|1 and the analytic continuation of eo for p —>■ p would 
coincide with the solution of Eq. H79|l with density p. As a consequence, we should 
have W > F{eo) = 1 in contradiction with the hypothesis W < 1. In conclusion, for 
pC^^min) = the characterization of the ground state in terms of W is the same as in 
the discrete case 

W > l4^eo< Wmin, (84) 
I^ < 1 ^ eo = «min. (85) 

Even when the potential levels V are continuous random variables, the definition 
H62|l and the relations H63|) and (|64|l remain unchanged, and by performing functional 
derivatives of the rescaled ground-state energy eo with respect to w, we have 

1 -1 



C{v) = 



Seo 



p{v) 



Sv {v — eo)^ 

Equation H86() holds for any density p{v] 
and we get 



p{v') 



dv' 



(V - eo)2 
Again, if p{v 



W <1, 
W>1, 



whereas for v > v„ 



C{v) 



p{v) 



:dv' 



(86) 

0, from Eqs. H84I85|I 
(87) 



W <1, 
W >1. 



(v - eo)2 [J {v' - eo) 

Finally, it is simple to extend Eq. (|68|l for the rescaled hopping energy in the 
ground state, 

-1 



hop 



p{v) 



:dv 



. • (89) 

Again, if p(wmin) = 0, wc see that whereas a non vanishing hopping energy is obtained 
for W > 1, the phase for < 1 is a frozen phase with eQ°^ = 0. 

We conclude this Section by discussing the results of numerical simulations on 
finite-size systems with a random potential having a continuous spectrum. The 
systems considered are spinless hard-core bosons with long range hopping in a one- 
dimensional lattice. Similar results, not shown here, arc obtained for spinless hard-core 
bosons and spinless fermions with first neighbor hopping. The relevant parameters 
are given in Table El No relevant differences are observed in two-dimensional lattices. 
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Figure 11. Universality of the tiiermodynamic limit of random potential 
systems: rescaled ground-state energy eo as a function of the rescaled potential 
strength 7/|i?Q'^'| for systems, Ixn.N, of A' spinless hard-core bosons with long 
range hopping in a one-dimensional lattice with n sites. The potential randomly 
associated with the states is y = ■yx, where x is a random variable in the interval 
[0, 1] with normalized constant density p{x) = 1. The solid line is the universal 
thermodynamic limit predicted by Eq. 1791 whereas the dots are the numerical 
results for the finite-size systems indicated in the legend. 

In Figs. II II and [T^ we illustrate the results for two distinct potentials. In both cases 
we have V = jx where a; is a random variable in the interval [0, 1] with normalized 
constant density p{x) = 1 in Fig. 1111 and with normalized linear density p{x) = 2x 
in Fig. Equation H79|l predicts, in the thermodynamic limit, a smooth rescaled 

ground-state energy cq = i?o/|£'o°''l ^'^^ potential in Fig. II II and a phase transition 
for the potential in Fig. El In the latter case, the singularity is localized by the 
condition W = 1, i.e. at 7/|i?Q°''| = 2. The results from the finite-size systems in Figs. 
Illl and lT^ clearlv show convergence towards the predicted thermodynamic behavior. 

Finally, in Fig. 1131 we show the rescaled hopping energy in the ground state, 
eg°'', as a function of the rescaled potential strength for the two potentials considered 
in Figs, nil and IT^ In the case of constant potential, we have a smooth transition 
between the values e^"^ = — 1 at 7 = and eg°^ — > for 7 3> |£^o''''|- On the other 
hand, in the case of the linear potential there is a discontinuity in the first derivative of 
eg°'' as a function of 7 at 7/|-Bo"''l = 2. This behavior is consistent with that observed 
in the simulations of finite-size systems. 

The random potential systems considered in this Section are a many-body 
generalization of the well-known Anderson models [H|, namely a single particle in a 
c?-dimensional lattice in the presence of a disorder potential. For the Anderson models 
with tight binding Hamiltonian the following rigorous results have been established 
1^. For d > 1, the eigenstates of the Hamiltonian undergo a localization phase 
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Figure 12. As in Fig. 1111 in the case of the potential V = -fX, where i: is a 
random variable in the interval [0, 1] with normalized linear density p(x) = 2x. 
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Figure 13. Rcscalcd hopping energy in the ground state ejj"^ as a function 

of the rescaled potential strength 7/|i?o^'| for the systems with constant and 
linear random potential considered in Figs. 1111 and 1121 The solid lines are the 
thermodynamic values predicted by Eq. 1891 whereas the dots are the numerical 
results for the finite-size systems indicated in the legend. 
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transition for large disorder; in one dimension, all states are localized for arbitrarily 
small non zero disorder. It is clear that our analysis does not apply to the single- 
particle Anderson models, in which the number of active links, i.e. the noninteracting 
ground-state energy Eq^\ remains constant in the limit M — > oo. Note that in the 
single-particle case the Fock dimension M coincides with the number of lattice points. 
However, in the many-body case the phase transition predicted by Eq. l(7??jl has some 
analogies with the phase transition in the Anderson models. The transition of the 
ground state to the frozen phase characterized by ejj"^ = takes place, universally, 
for < 1. In [Appendix B| we show that this condition on the functional 
is equivalent to a condition of large disorder, i.e. large support of p{v) or large 
potential strength. Moreover, as predicted by Eqs. lP?7| and (p^ and also checked 
in the numerical simulations for finite-size systems, in the frozen phase the ground 
state of the system \Eo) condensates into a superposition of the Fock states which are 
eigenstates of the potential operator with eigenvalue V ~ Knin- This condensation, 
which is a sort of superlocalization in the Fock space, does not correspond, in general, 
to a localization in the lattice space and takes place for any dimension d. 

8. Conclusions 

We have revisited the probabilistic approach [SJ (HJ [7j to evaluate the ground state of 
general matrix Hamiltonian models from a different point of view. The asymptotic 
probability density of the potential and hopping multiplicities which is the core 
of the approach is considered known and an exact equation for the ground-state 
energy Eq is obtained in a closed form. For a class of systems, which includes 
the uniformly fully connected models and the random potential systems, we have 
that the above mentioned probability density factorizcs into separated potential and 
hopping contributions, the potential part being of multinomial form. This permits to 
write a very simple scalar equation which relates Eq to the lowest eigenvalue of the 
Hamiltonian operator with zero potential iSg"'' . Such relation becomes exact in the 
thermodynamic limit in which iJg"'' diverges. 

The factorization of the probability density into separated potential and hopping 
contributions has a different origin in the uniformly fully connected models and in the 
random potential systems. In the former case, due to the particular structure of the 
hopping operator which connects with equal probability any state to any other one, 
the potential values assumed by the system during its evolution are uncorrelated to 
the number of allowed transitions, equal for any state. In the latter case, the nature 
of the potential operator ensures that the potential levels are associated randomly 
with the states of the system. Correlations between potential and hopping values are, 
however, reintroduced by the dynamics at a degree which depends on the structure of 
the hopping operator. The correlations disappear if the number of allowed transitions 
from any state diverges, i.e. Aav ~ l^'o"^! ~* This explains why the relation that 
we have found between Eq and Eq^^ becomes exact in the thermodynamic limit. 

In the thermodynamic limit, we find that the rescaled energy eg = Eq/IEq^'^I is a 
universal function of the potential levels, rescaled by |£^o°''l' ^^'^ corresponding 
degeneracies, for the uniformly fully connected models, or probabilities, for the random 
potential systems. In the case of the random potential systems this means that eo 
does not depend on the nature of the hopping operator. In general, if the degeneracy 
(probability) associated with the lowest potential level vanishes, the ground state of 
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the system undergoes a quantum phase transition between a normal phase and a 
frozen phase characterized by zero hopping energy. The control parameter W of the 
phase transition is related to the whole set of levels and corresponding degeneracies 
(probabilities) of the potential. This parameter measures the inverse amount of the 
spread (disorder) introduced by the potential operator (random potential). In the 
frozen phase, corresponding to large spread (disorder), i.e. small W, the ground state 
condensates into the subspacc spanned by the states of the system associated with the 
lowest potential level. 

We have performed numerical simulations on different finite-size random potential 
systems. In particular, we considered many-body lattice systems in which the particles 
are fermions or hard-core bosons, with long range or first neighbor hopping, in one or 
two dimensions. The results found for the ground state are nicely consistent with the 
universal behavior predicted in the thermodynamic limit. 
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Appendix A. Proof of Equation (|32|) 

Let us rewrite the probability density VNiNu) in terms of its Fourier transform VNio) 
Pn{Nu) ^ I dq e^o^.[■PN(q)]~^N(q,u)^ -^^ 



If we indicate with {v^-^^ . . . with ai, . . . , £ , a generic cumulant of order 

fc, we have 

log7'A.(g)=log(e'^(-'^))^ = 5:^((z.,zq)'=)^ . (A.2) 

fe=i 

For any given TV, due to the inequalities {^ai ■ ■ ■ fJ-Uk) ^ — valid for any fc, and 
due to the asymmetry 'Piv(/j) ^ 7'Ar(— //), the series in Eq. (|A.2p converges for every 
q G (see, for example [TO]^. 

Let us introduce the rescaled cumulants of order fc, in a compact notation S^^'''), 
which are defined as the tensors of rank fc with components 

- ('^"i ■.■vc.k)%\ ai, . . . , «fc e J^, (A.3) 

(A;) 

and let I]ai....,afc be their asymptotic values 

Si';l..,Q, - Hm 7V^-i {v^, . . . v^,)%^ ^ hm 1 (iV„, . . . . (A.4) 

These limits exist and are finite since the irreducible and finite Markov chain formed 
by the evolving configurations has a finite correlation length Nc with respect to the 
number of jumps N. As a consequence, up to corrections exponentially small in N/Nc, 
for N ^ oo we have 

fc=l ■ aie.y^ akGJ^ 

= Ng{q), (A.5) 
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where the function g{q) is independent of N. 

By using the result (|A.5() and evaluating the integral in Eq. l)A.l|l by steepest 
descent, we find the asymptotic logarithm equality 

Pn{Nu) ~ exp [Nlj{u)] , (A.6) 

where ijj{v) — giq^^) — i{q^^ , i^) is the function at the exponent evaluated in the saddle 
point q^P(i/). Equation (|32|1 follows immediately. 

Appendix B. The critical condition W = 1 in terms of the disorder 

In this Appendix wc show that the condition under which the ground state of the 
system undergoes a quantum phase transition, namely = 1, is equivalent to 
WA = 0(1), where A is the amount of disorder defined as A — v — fmin- 

Let us consider a density p(v) with p(w,„i,i) ~ and suppose that all the statistical 
moments of p{v) are finite, (u*^) < cxo, with k non negative integer. Let us also 
suppose, for simplicity, that the density p{v) is symmetric around the mean value v so 
that all the odd centered moments are zero and, furthermore, for any v in the support 
of p{v) we have 

\v -V\ < \V-V„un\- (B.l) 

By using the geometric series, which is convergent due to the above strict inequality, 
for any small positive S we have 

-dv^ I /\' dv+- — / L^dv.{B.2) 



^mi n Jv-i-n^,^ ^ ni i II ^ ^ m i ii 



1-5 ^,^0 ~ "^'min) 



In the interval [wmin + S, Wmax] , for sufficiently small positive S wc also have 

_ = 1-- < 1. (B.3) 

Equation IjB.Sjl implies that the series in the following inequality, 

°° /■"max f —\2k OQ / ^ \1k 

E/ y ~ dv < - ^S)Y^il- , (B.4) 

converges absolutely for any positive (sufficiently small) S, so that in the second term 
of the r.h.s. of Eq. ljB.2|l we can exchange the order in which the operations of 
integration and series appear. Finally, by performing the limit 6^0 and using the 
fact that p(umin) = 0, we get 

W= ._ ^ L^f. (B.5) 

(V Wmin) {v - Umin) 

From Eq. IjB.Sp we see that W is a. decreasing function oi A ^ v — I'nun and 
W = 0{1)/A, so that the critical condition amounts to A = 0(1). By recalling 
that V = V/IEq^'^I, in terms of the non rescaled potential the critical condition reads 
V — Knin — 0(|i?Q°''|). This result is a generalization to many-body systems (in 
any dimension) of the result established for the single-particle Anderson models with 
dimension d> 1 |S1 ■ 
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